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Abstract. The essential spectrum of multiparticle Brown- 
Ravenhall operators is characterized in terms of two-cluster de- 
compositions for a wide class of external fields and interparticle 
interactions and for the systems with prescribed symmetries. 
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1 Introduction 

It is well known that the eigenvalues of the one-particle Dirac operator are 
in much better accordance with the spectroscopic data then the eigenvalues of 
the Schrodinger operator. However, due to the presence of the negative con- 
tinuum of positronic states the many-particle Coulomb-Dirac operator has no 
eigenvalues and its essential spectrum is the whole real line. Coupling with the 
quantized electromagnetic field does not correct this situation. However, there 
are ways to construct a semibounded operator which will take the relativistic 
effects into account. Such models, although nonlocal, find their applications 
in numerical studies of heavy elements and cosmology, where the relativistic 
effects cannot be ignored. 

The most obvious choice of the kinetic energy (sometimes called Chandrasekhar 
or Herbst operator) given by yp^^ + m^c?, p and m being the momentum 
and mass of the particle, suffers from the lack of semiboundedness for nuclear 
charges exceeding 87, as shown in [8]. Most other operators considered in the 
literature are obtained by reducing the (multiparticle) Dirac operator onto some 
subspace on which it becomes semibounded. One of such models, extensively 
studied recently, is by Brown and Ravenhall [1], see also Bethe and Salpeter 
f3] , Sucher [Tni [T7] . In this model we require that every particle stays in the 
positive spectral subspace of the free Dirac operator. Since the multiplication 
by interaction potentials does not leave this subspace invariant, the potential 
energy terms should be projected back by the corresponding projector. 
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The mathematical study of the Brown-Ravenhall operator started from the 
one-particle case in the article of Evans, Perry, and Siedentop 6J. The au- 
thors have proved that the atomic Hamiltonian is semibounded from below for 
nuclear charges not exceeding 124. This makes the Brown-Ravenhall model 
applicable to all existing elements. It was also proved in [B] that the essential 
spectrum of the one-particle atomic Brown-Ravenhall operator is [mc^, oo) 
with m being the mass of the particle, and that the singular continuous spec- 
trum is empty. 

Further studies of the Brown-Ravenhall operator include the improved lower 
bounds by Tix [191 ISO] (see also Burenkov and Evans [5] ) in the atomic case, 
the proof that the eigenvalues of Brown-Ravenhall operator are strictly bigger 
than those of the one-particle Dirac operator by Griesemer et al. [7], proofs 
of stability of one-electron molecule by Balinsky and Evans [2] and the proof 
of stability of matter by Hoever and Siedentop |2] . The essential spectrum of 
the many-particle operator was characterized by Jakubai3a- Amundsen [lOllllj . 
and Morozov and Vugalter [14] in terms of two-cluster decompositions. This 
is usually referred to as HVZ theorem after the well known result for the many 
particle Schrodinger oprator. It is also shown in [T3j that the neutral atoms or 
positively charged atomic ions have infinitely many bound states. 
In all these previous studies the nuclei were considered as fixed sources of the 
external field, the particles were assumed to be identical, and the interaction 
potentials were purely Coulombic. 

In this paper we generalize the HVZ theorem of [TOlUlllTl] as follows: We allow 
any number of (massive) particles of the system to be identical. We allow quite 
general matrix interaction potentials. In particular, our result applies in the 
presence of the magnetic fields if the vector potential decays at infinity in some 
weak sense. Another problem we address is the reduction to any irreducible 
representations of the groups of rotation-reflection symmetry and permutations 
of identical particles. Note that such a reduction allows to analyze the eigen- 
values of some irreducible representations even if they are embedded into the 
continuous spectrum of another representation. Existence of such embedded 
eigenvalues is well known for atomic and molecular Schrodinger operators. 
From the technical point of view, the nonlocality of the model due to the 
presence of the spectral projections of the free Dirac operator is overcome with 
the same ideas as in [14]. One more complication should be stressed: for the 
Brown-Ravenhall operator the center of mass motion cannot be separated in 
the same way as it is usually done for Schrodinger operators, where the complete 
Hamiltonian without external field can be represented in suitable coordinates 
as 

n^A^i + i(g)B, 

where A describes the free motion of the center of mass and B is the internal 
Hamiltonian of the system (see [H]). Such a decomposition appears to be 
especially fruitful in the presence of rotation symmetries. Since it cannot be 
obtained for pseudorelativistic operators due to the form of kinetic energy. 
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we have used completely different approach based on the commutation of the 
Hamiltonian with the absolute value of the total momentum of the system. 
Note that the proof of the HVZ theorem for a system of particles described by 
the Chandrasekhar operator, was till now not known in presence of rotation- 
reflection symmetries (see the article of Lewis, Siedentop and Vugalter [T3] 
for the proof without symmetries). Such a proof can now be obtained as a 
simplified modification of the proof given in this paper. 

In Section [5] we introduce the model and make the necessary assumptions. At 
the end of this section we formulate the main result in Theorem [6l The rest of 
the article contains the proof of this theorem. 



2 Setup and Main Result 



[A, B] = AB — BA is the commutator of two operators. (•, •) and || • || stand for 
the inner product and the norm in L2{M.^'^,C'^''), where d is the dimension of 
the underlying configuration space. Irrelevant constants are denoted by C. Iq 
is the indicator function of the set fi. For a selfadjoint operator A we denote its 
spectrum and the corresponding sesquilinear form by cr(A) and {A-, •) ~ (•, A-), 
respectively. We use the conventional units h = c = 1. Sometimes we denote 
the unitary Fourier transform by ^. 

In the Hilbert space L2{M.^^, C"*) the Dirac operator describing a particle of mass 
m > is given by 

Dm ~ —ioL ■ V + /3m, 

where a := (ai,Q;2,Q;3) and (3 are the 4x4 Dirac matrices [18]. The 
form domain of D„i is the Sobolev space i/i/2(M3,C4) and the spect rum IS 
(— (X), — m] U [to, +cxd). Let A„i be the projector onto the positive spectral sub- 
space of Dm'- 

" - 2 + 2V-A + TO2 • ^'-'^ 

We consider a finite system of N particles with positive masses to„, n — 
1, . . . ,iV. To simplify the notation we write Z?„ and A„ for -Dm„ and A^^, 

N 

respectively. Let fyjq := ® A„L2(IR."^, C**) be the Hilbert space with the inner 

n— 1 

product induced by those of ® L2(IR^ C^) ^ L2(M3^, C^"). In this space the 

n—l 

A''-particle Brown-Ravenhall operator is formally given by 

, N N s 

nN^A''[J2iDn + V,,) + J2Un,)A'', (2.2) 

with 



n<j 



^ N 



A^ IT ^" = ® (2-3) 



n-1 
n—l 
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Here and below the indices n and j indicate the particle, on whose coordinates 
the corresponding operator acts. In (|2.2p Vn is the external field potential 
for the n"* particle, i.e., the operator of multiplication by a hermitian 4x4 
matrix-function V^(x„), n = 1,...,N, and Unj is the potential energy of 
the interaction between the rt"* and j*'' particles, given by the operator of 
multiplication by a hermitian 16 x 16 matrix-function [/„j(x„ — Xj), n < j = 
1, . . . ,N. More explicitly, if we let sj S {1, 2, 3, 4} be the spinor index of the 
j*'* particle, then 

(14V')(xi, si; . . . ; x„, s„; . . . ; x^v, sn) 
:= ^T4'*"''*"(x„)?A(xi,si; . . . ;x„,s„; . . .■,xn,sn), 

and 

: Sn ; • • • ; Xj , ; . . . ; Xjv , S Af j 

:= ^ C/^^^"'*"'''(x„ - Xj)'0(xi,si; . . . ;x„, s„; . . . ; x^, Sj-; • • ■ ;xAr, sw). 

Before we make other assumptions on the interaction potentials, let us consider 
possible decompositions of the system into two clusters. Let Z — {Zi, Z2) be a 
decomposition of the index set / :— {1, . . . , N} into two disjoint subsets: 

I = ZiU Z2, Zi n ^2 = 0- 

Let 

N,:=#Z„ J = 1,2 (2.4) 

be the number of particles in each cluster. We will write n^ffj if n and j belong 
to different clusters. Let 

nz,i + "^-) + E ^"J-' (2.5) 

n<j 

nz,2 ■.= Y,Dr,+ (2-6) 

n<j 

We omit Ti-zj Zj = 0, j = 1,2. Let us introduce the operators corresponding 
to noninteracting clusters, with the second cluster transferred far away from 
the sources of the external field: 

nz.r-^^z.jnz.jAzj, m f)z,r.= ® A„L2(K^C4), j = l,2, (2.7) 
where 

neZj 

We make the following assumptions: 
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Assumption 1 There exists C > such that for any Z and j — 1,2 
|(Hz.j(^»| C||(p||jji/2||^A||hi/2, for any (p, ^ e (g) H^^^{R^,C^). (2.8) 

neZj 

For Coulomb interaction potentials p.8|) follows from Kato's inequality. 
Assumption 2 There exist Ci > and C2 G K such that for any Z 



neZj 

for any € ® A„i^l/2(M^C'^), j = l,2. 
Remark 3 Note that for ^ e (g) A„i?i/2(R3, C"') the metric 

n GZj n G Zj 

is equivalent to the norm of ijj in (g) H^^^(R^,C^), since 

riGZi 



(2.9) 



A„AiA„ = A„|Ai|A„ = A„v/-A + ?7i2A„. (2.10) 

An equivalent formulation of Assumption [2] is that the operator Ti.z,j is semi- 
bounded from below even if we multiply all the interaction potentials by 1 + e 
with e > small^enough. This is only slightly more restrictive than the semi- 
boundedness of 7iz,j ■ 

Assumption 4 For any R > there exists a finite constant Cr ^ such that 

^ / r \ 1/2 w / /. \ 1/2 

V / |K(x)|'dx +V / |C/„,(x)|'dx i^CB. (2.11) 

This means that the interaction potentials are locally square integrable. 

Assumption 5 For any e > there exists R > big enough such that for all 
n<j = l,...,N 

l|K/{|x„|>fl}V^IK£|||£'n|'/'^||, for all ^^H^'^W^X"), (2.12) 

and 

l|f/«,/{|x„-x,|>i?}^|| <emin{|||i5„|i/2^||,|||i?,f/V||}, 13) 

for all e i^^/2(M^Cl6). 
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By Remark [3] this assumption is weaker then the decay of Loo norms of the 
interaction potentials at infinity. 

It foUows from (|2.9p and Remark[3]that for any Z there exists a constant C > 
such that for any V' € <S) A„i/i/2(R'\ C"*) 

n<£Zj 

+ 11^11'), J = 1,2. (2.14) 

Hence by Assumptions [1] and [21 the quadratic forms of operators (|2.7p 
(and, in particular, TIn) are semibounded from below and closed on 
(g) A„i/i/2(M3^C'*). Thus these operators are well-defined in the form sense. 

Some particles of the system (say, k*'^ and /*'') can be identical (in which case 
iTT-k = mi, Vk = Vi, and Ukj = Uij for all j). Then the operator Hn can be 
reduced to the subspace of functions which transform in a certain way under 
permutations of identical particles. The most physically motivated assumption 
is that any transposition of two identical particles should change the sign of 
the wave function tp E S)n describing the system. This is the Pauli principle 
applied to the identical fermions (the model describes spin 1/2 particles, thus 
fermions). 

Let n be the subgroup of the symmetric group Sn generated by transpositions 
of identical particles. We denote the number of elements of 11 by hu- Let E be 
some irreducible representation of H with dimension cIe and character ^g. For 
V' G iOAT let 

pE^,^^Y^J^^^^ (2.15) 

where tt is the operator of permutation: 

(7rV')(xi, si; . . . ; xat, sn) = ■(/'(x7r-i(i), s^-i(i); . . . ; x^-i(jv), s^-i(Af)). 

Here si, . . . ,sn are the spinor coordinates of the particles. The operator 
defined in (|2.15p is the projector to the subspace of functions in which trans- 
form according to the representation E of H. Since any tt g H commutes 
with Hn, reduces Ti-N- Let Hf] be the corresponding reduced selfadjoint 
operator in 

For a decomposition Z = {Zi, Z2) let H^ be the group generated by transposi- 
tions of identical particles inside Zj, j = 1,2. For any irreducible representation 
Ej of Hj with dimension and character the projection to the space of 
fimctions in Sjj transforming according to Ej under action of H? is given by 



Essential Spectrum of Brown-Ravenhall Operators 7 



where /ijjz is the cardinaUty of IlJ. Projectors P^^ reduce operators H-zj- We 



introduce the reduced operators in 



:=P^^i3z,, , J = 1,2. 

Given an irreducible representation _B of 11 and a decomposition Z = {Zi, Z2), 
we have 

c e {S^z\^^t2)^ (2-16) 

[El ,£2) 

where Ei^2 are some irreducible representations of IlJ. We write {Ei, E2) -< E 

if the corresponding term cannot be omitted on the r. h. s. of (|2.16p without 
violation of the inclusion. 

Apart from permutations of identical particles the operator Ti^ can have some 
rotation-reflection symmetries. Let 7 be an orthogonal transform in R^: the 
rotation around the axis directed along a unit vector through an angle ip^ , 
possibly combined with the reflection x l—^ — x. The corresponding unitary 
operator Oj acts on the functions V' G as (see [TB], Chapter 2) 

Af 



Here S„ = — A is the spin operator acting on the spinor coordinates of 
the n*'' particle. The compact group of orthogonal transformations 7 such that 
commutes with Vn and Unj for all n, j = 1, . . . , TV (and thus with Ti.^) we 
denote by T. Further, we decompose S^ff into the orthogonal sum 

^N^ (B iD^"-'', (2.17) 

where Sj'^" '^ consists of functions which transform under Oj according to some 
irreducible representation Da of T, and A is the set indexing all such irre- 
ducible representations. The decomposition (|2.17p reduces Ti-f^. We denote the 
selfadjoint restrictions of to ^^"'^ by H^"'^ . For any fixed irreducible 
representation D with dimension do and character the orthogonal projector 
in fjN onto the subspace of functions which transform according to D is 



1'- 
r 



where fi is the invariant probability measure on F. 

For j — 1,2 let Dj be some irreducible representations of F with dimensions 
and characters Coj- The corresponding projectors in ^)z,j are given by 



r 
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where O^j is the restriction of 0~f to S)z.j'- 

Given representations Dj and Ej, projector P^o pEj — pEj pDj j-educes Hzj- 
We denote the reduced operators in 

^zf' ■■^P'''P'''^z,j 

hyH'^f\Let 

^,iZ,D„E,) ■.= Mainzf'). (2.18) 
We write {Di, Ei; D2, E2) -< {D, E) if the corresponding term cannot be omit- 
ted on the r. h. s. of 

(Di ,Ei ) 

without violation of the inclusion. For Z2 7^ let 
>i{Z,D,E) 

^(mi{>ci{Z,Di,Ei) + >C2{Z,D2,E2) : {Di, Ei; D2, E2) < {D,E)}, Zi ^ 0, 

'~\>t2{Z,D,E), Zi-0. 

(2.19) 

The main result of the article is 

Theorem 6 For N Cz N let D be some irreducible representation of T, and E 
some irreducible representation ofH, such that p^p^ ^ Q. Then 

CTessO-i^'^) = [>c{D,E), 00), 

where 

x{D,E)=inm{^(Z,D,E) : Z = (^1,^2), Z2 + 0}. (2.20) 

Remark 7 We only need Assumption\^for the operators Ji^'-^' which appear 
m (j2T8)l . (j2T9)l . 

3 Commutator Estimates 

3.1 One Particle Commutator Estimate 

Lemma 8 Let x G (^^(R-^) (i. e. a twice- differentiable function with bounded 
derivatives) . Then for nin > the norm of the operator 

[X,A„] :L2(M^C4)^i^l(M^C4) 
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satisfies 

^C{m^){\\^x\\L^ + \\d\U^). (3.1) 

Here\\d^x\\L^= max |a|;x(z)|- 
fcJe{i,2,3} 

Proof. In the coordinate representation for / £ Cq(R'^,C'*) the operator A„ 
acts as 

(A„/)(x).M + |^^H,^^ / ^^^^^iKlx-y|)/(y)rfy 

y-x|>e 

+ 4;^y(,^ |x-y| + |x-yp i^oK|x-y|)j/(y)dy, 

where the hmit on the r. h. s. is the hmit in L2{M.^, C**) (see Appendix B of [l4] . 
where this formula is derived in the case to„ = !)■ The rest of the proof is an 
obvious modification of the proof of Lemma 1 of [14 , where the case m„ = 1 
is considered. • 



Remark 9 Since we only deal with a finite number of particles with positive 
masses, we will not trace the m-dependence of the constant in ()3.ip any longer. 



3.2 Many Particle Commutator Estimate 

Lemma 10 For any d,k £ N there exists C > such that for any % G 
anduG H^^^{R'^,C'') 

llx"llHi/2(Rd^c'=) ^ C'(llxllLoo(a'') + II^xIIloo(R''))II"IIhi/2(R'',c'=)- (3-2) 

Proof of Lemma [TUl We choose the norm in iJ^/^(M'', C*^) as (see [1 , The- 
orem 7.48). 

II ii2 _ II ii2 _L f f l"W ""(y)l^ ^ . 

ll"ll_Hi/2(Rd,c'=) II'"IIl2(R'',C'=) + J J jjj. _ y\d+l "^"y- 

Then 

|2 



II ii2 II 112 ^11 x(x)m(x) - x(y)M(y) ^ , 



^11 ii2 II 112 , / / / kW| hW-w(y)| , |x(x)-x(y)| \u{y)\ , 

^\\x\\lju\\t + II (J \^-y\d+l + |x„y|d+l ^l^'^'^y 



^ II l|2 II ||2 _L / |X(X) -X(y)r ^ II ||2 

< WxWLjmlm/^ + sup / \x-y\d+i "-^Ml-- 



(3.3) 
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The supremum on the r. h. s. of (|3.3p can be estimated as 



/■|xW-x(y)| , ^ f |xW-x(y)| , 

sup, / ix-vK+i '^'^ ^ '""P, /, , Ix-vK+i '^'^ 

ygRdJ |X y| yGR'' J|x-y|^l \^ Yl 

+ SUP / kW^^dx<|§'^-i|(||Vx||L+4||x||L), 

yeK'' J|x-y|>l |X — y| 



where is the area of {d— 1) -dimensional unit sphere. Substituting (13. 4p 

into (13.31) we obtain (13.21). • 



Lemma 11 For any x G C^(M'^^) the operator [%, A^] zs bounded in 
i7i/2(M3^,C4"), and /or any V e -^^/^(K^^, C*" ) we have 

(3.5) 

with C depending only on N and the masses of the particles. 

Proof. Successively commuting x with A„, n = 1, . . . , TV (see (12. 3p ) we obtain 

N n-1 N 

[^'A^]=EnAfe[x,A„] n (3.6) 

n— 1 k—1 l—n+1 



where the empty products should be replaced by identity operators. By (|2.ip 
the operators A„ are bounded in H^^^ for any n — 1, . . . ,N. This, together 
with p.6p . and Lemmata [51 and [TUl implies (|3.5p . • 



4 Lower Bound of the Essential Spectrum 
In this section we prove that 



Ma,,s{n^^'')> x{D,E). (4.1) 
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4.1 Partition of Unity 

Lemma 12 There exists a set of nonnegative functions {xz} indexed by possi- 
ble 2-cluster decompositions Z = [Zx^Z^] satisfying 

1. Xz ^C"°{^'^'^) for all Z; 

2. xz(«X) = X2(X) for all |X| = 1, > 1, Z2 ^ 0; 

3. 5]xl(X) = 1, for all X e M^^; (4.2) 

z 

There exists C > such that for any X G supp yz 

4. . , (4.3) 
mmjlxj - x„| : G Zi, x„ e Z2; |x„| : x„ G ^2} > C|X|; 

(7x1 , ■ ■ • , ) = Xz (xi , . . . , xat ) /or any orthogonal 

5. 

transformation 7; 
6. Xz is invariant under permutations of variables preserving Zi,2- 

Proof. The proof is essentially based on the modification of the argument 
given in [TS], Lemma 2.4. 



1. We first prove that for any X = (xi, . . . ,xn) G with |X| = 1 there 
exists a 2-cluster decomposition Z = {Zi, Z2) such that 

min {|xj - x„| : x^ G Zi, x„ G Z2; |x„| : x„ G Z2} > iV"^/^ (4.4) 

Indeed, let k be such that |xfc| ^ \ for all j = 1, . . . , TV. Then, since |X| = 1, 

|xfc|^iV-i (4.5) 

Choose Cartesian coordinates in R'^ with the first axis passing through the 
origin and x^, so that x^ = (|xfe|,0, O). Consider N regions 

Ri := {xgK^ M/N}, 

Ri ■■= {xgM^ e ((/-l)|x,|/iV, /|xfe|/7V]}, l = 2,...,N. 

At least one of these regions does not contain Xj with j ^ k. Let be the 
maximal index of such regions. Let Z2 be the set of indices n such that x„ G 
U Ri. Z2 is nonempty since x/c G Z2. Setting Zi / \ Z2 we observe that 

/>/o 

min {|xj - x„| : Xj G Zi, x„ G Z2; |x„| : x„ G Z2} > |xfe|/iV, 
which together with (|4.5p implies (|4.4p . 



2. Choose r] G C°°(R+) so that 

r,{t) = 



0, tG[0,l] 

1, tG[2,oo). 
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Let 

«x,Jn;<l^)n''(p^)"(^M)^ 

I n£Z2 

[l-77(2|X|), Z2 = 0. 

(4.6) 

Functions satisfy conditions 1, 2, 4 (with C = iV"^/^), 5, and 6 of 

Lemma [121 Moreover, by the first part of the proof 



^Cz(X)^l, foraU Xe 



p3N 
Z 



Hence all the conditions are satisfied by the functions 

z 



xz -cy^(Ec.) 



Let 

xf(X) :=Xi?(X/i?), (4.7) 

where the functions xz are defined in Lemma [T21 The derivatives of decay 
as i? tends to infinity: 

livxf Hoc ll^^xf llo, Ci?-2. (4.8) 

To simplify the notation we omit the superscript i? further on. 
4.2 Cluster Decomposition and Lower Bound 

We now estimate from below the quadratic form of Ti.^'^ on a function tfj from 
^'^'^ n i i/i/2(M^ C^), which is the form domain of H'^'^ . 

. N N s 

(<'^^, ^) = ( E(A. + Vn) + E Un, ] E xlV', ^) 
^n=l n<j ^ Z 

. N N s 

= EU E(^" + + E ) XZ^, XZ^>- 

Z ^n=l n<j ' 

Here we have used (|4.2p and the relation 

JV JV N 2 

E(/' E ^nixlg)) = E(x^/' E ^nixzg)) + E(/' E (^)5> (4-9) 



Z n=l Z n=l Z 11=1 
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N 

which holds for any f,g £ » i7^/^(M^, C*). The last term on the r. h. s. of 

n— 1 

is equal to zero due to (14. 2p . Thus 

Z=(Zi,Z2) ^ 

+ {{Hz.i + Hz,2)[xz, A^]V, A^'xzi^) 



n£Z2 n<j 

The terms at the last line of (|4.10p can be estimated as 

(E v;xl^,V') + (E/^n,xlV',V'> + II V^ir) (4.11) 

nGZ2 n<j 

with ei{R) ^ as i? -> oo, due to ^J^, ^J^, g3D, g21), and ((m)) . The 
terms at the second and third lines of (|4.10p can also be estimated as 

{{nz.i + Hz^2)[xz, A^]^, A^xzV^) + {{nz.i + nz.2)xzi^, [xz, A^]V) 
^ -e2(i?)((<'''^,V^) + IIV^IP), e2{R) 0, 

due to ([^ . ([X^ . and In the case Z2 7^ we estimate the 

terms at the first line of (|4.10[) in the following way (recall the definitions of 
>Cj{Z,Dj,Ej) and >c{D,E) given in (PTS)) and ^^TIU^i ): 

((Hz.i+W2,2)A^Xz^,A^Xz^) 

J2 ((Wz,!^''^^''^ +^z,2P''^P^^)A^XzV',A^XzV') 

{Di,Ei;D2,E2)-<{D,E) 
(Di,Ei;D2,E2)-<{D,E) 

+ K2{Z, D2, £;2)P^^P^^) A^xzV-, A^xzV-) 
^ >c{D,E){A^XziJ,f^''xz^) 
= ><(I?, i?)(x|^, ^) + M^, i?)([A^, xz]'^, XzV') 
+ >^(Ai?)(A^XzV',[A^,Xz]V')- 

(4.12) 

By (|3^ . (|33)) . (|48l) . and ((2T4| the last two terms on the r. h. s. of (|4T2)) can 

be estimated as 

>c{D, i?)([A^, xz]^, XzV-) + >^(^, i?)(A^xzV', [A^, Xz]^) 

^-e3(i?)((<'''^,V^) + 11^-11'), e3(i?)--0. ^^-^^^ 
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Substituting the estimates (|4.1ip — (|4.13p into (|4.10p we obtain 

-s4R){{n'^'''i^,i^) + ei{R) — > 0. 

(4.14) 

4.3 Estimate Inside of the Compact Region 

It remains to estimate from below the quadratic form of the operator TL^'^ 
on the function A^X(/,0)'0- Note that according to Lemma [T^] and (|4.7p 
suppx(7,0) C [—2R,2R\^'^. To simphfy the notation let 

Xo := X{i,0)- 

Lemma 13 For M > let 

Wm ■■= {p e : \p^\ ==; Af, 1 = 1, . . . , 37V}, Wm \ Wm- 

There exists a finite set Qm ^ L2{M-^^) such that for any f G L2{M.^^) with 
supp/ C [-2i?,2i?]3^, fJ-QM holds 

The proof of Lemma [12] is anafogous to the proof of Theorem 7 of [21] and is 
given in Appendix C of [14j . 

It foUows from ([^ that for any M > 

N 
n=l 

(4.15) 

Here -fpf?^^^ is the operator of multiphcation by the characteristic function of 

Wm in momentum space. 
We choose 

M ■.= 8{>c{D,E) + C2)Cj;'^ (4.16) 

and assume henceforth that / :— xoV-' is orthogonal to the set Qm defined 
in Lemma 1131 Since in momentum space the operator Dn acts on functions 
from A„L2(IR'^,C) as multiphcation by a/|pP + , by construction of Wm 
we have 

N 

(^Z?„/^^^A^Xo^,A^Xo^) > A/||/^^A^xo^f. (4.17) 

ra=l 



Essential Spectrum of Brown-Ravenhall Operators 



15 



Inequalities (|4.15p and (|4.17p imply 

^CiM(|1/^^^XoV'II - ||/^^,^[A^,Xo]V'||)'-C2||xo^||' 

- S^i^, i?) + C2) II [A^, XolV-f - CallxoV-f • 



At the last step we have used ()4.16p . The second term on the r. h. s. of ()4.18p 
can be estimated analogously to (14. 13^ as 

-%{>i{D, E) + C2) II [A^, xol^f ^ -e5(i?)((H^'''V', V-) + IIV-lP), £5{R) 0. 
For the first term on the r. h. s. of (|4.18p Lemma [T51 implies 

4|l^w,,Xo^|P > llxoV-lP. (4.19) 
As a consequence of (|4.18l) — (|4.19p . we have 

(<^^A^XoV',A'^XoV'>^MAi?)llxoV'li'-e5(i?)((<'''^,^) + ||Vll'), 



£5(i?) 0. 



(4.20) 



4.4 Completion of the Proof 
By HUM . Km . and (OD 

>M^,i?)IIV^II'-£6(i?)((<'''^,V') + 11^11'), ee{R) ~- 0. 

for any in the form domain of Ti^'^ orthogonal to the finite set of functions 
(cardinality of this set depends on R). This implies the discreteness of the 
spectrum of 7i^' below >c{D, E) and thus (|4.ip . 

5 Spectrum of the Free Cluster 

In this section we characterize the spectrum of the cluster Z2 which does not 
interact with the external field. 

Proposition 14 For any irreducible representations £'2, £-2 0/ rotation- 
reflection and permutation groups the spectrum o/7i^2'^^ is 

with some >C2{Z, D2, E2) G K. 
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Proof. Let us introduce the new coordinates in the configuration space . 
of the ciuster Zi = {ni, . . . , nj^^^ in the same manner as it is done in [13] . Let 
M := J2neZ2 totai mass of the particles constituting the cluster. 

Wc introduce 



n". (5.1) 



■Tlk + l ^11 1 J 



The Jacobian of this variable change is one. yo is the coordinate of the center 
of mass, whereby y^, k = 1, . . . , iV2 — 1 are the internal coordinates of the 
cluster. Accordingly, 



x„i = yo - rrin^+.yk, 
fe=i 
^ N2-1 

^n,+i = yo + y/ - T7 "^"f^+iyfe' ^ = I,- - - ,^2 - 1- 



(5.2) 



M 
fc=i 



The momentum operators in the new coordinates are 

N2~l 

^'^ t[ (5.3) 
Pn, :--*Vx,.^ =^P + (-*Vy,_J, fc = 2,...,7V2, 

where P is the total momentum of the cluster: 

neZ2 

Let J-Q be the partial Forurier transform on 2'^^ defined by 

(J^o/)(P,yi, ■ ■ ■ ,yjV2-i) := I /(yo,yi,---,yiV2-i)e~*^^Vyo. 



By p.6p we have 
where in the new coordinates 



^2-1 

''^2,2^' — X! +''^""*")+ X! X! t^nfcn,(yfe-yi): 

nGZ2 k=2 l<fe</^Af2-l 

(5.4) 

Az,2 := n (5.5) 

nGZ2 
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p„ are givn by (|5.3p , and P should now be interpreted as multiplication by the 
vector-function. The operators (|5.4p and (|5.5p obviously commute with *P := 
|P|. The operator J^q^^J^q (unlike JFq^^PJFq) is well-defined in ij^j' ^' since 
it commutes with P^^ an in i3z^2- This implies that Ti.^^^^ commutes 
with Tq'^^To- 

Let uj P/*P G S*^. We decompose the Hilbert space S)^^^^^ into the direct 
integral 

^z'f' = / ®S,^r'''^rd^. (5.6) 

The fibre space Sj^^^^'^ '^ can be considered as a subspace of L2{M.^^'^~^ x 
S'^jC'*"^) with the inner product 



JK3("2-1)xS2 

For / e ^z^2 ^ corresponding element of 2' ^ is given by 

/ip := ^o/||p|=?p- 

We have 

/•OO 

11/11'- / ll/'pll**P'c«'P (5.7) 
Jo 

in compliance with (|5.6p . The form domain of 7i^2'^^'^ is 

:= A|2F^^P^^ffi/2(M3(^^-i) X C^"^), 

where 2 is given by (|5.5p with the only difference that we should replace P 
by uj^ in (|5.3p . The operators on fibres of the direct integral (|5.6p are 

where 'H^'^.f''^'^ is given by the r. h. s. of (|5.4p with P replaced by a;*P in (|5.3p . 
We thus have 

/•OO 

7^22'^' = / ®H^^2'^'''^?52d<p. (5.8) 



The spectrum of 7i^2'^^ can be represented as 



a(H^-''^)=ess U c7(H^-^-'P), (5.9) 
<peR+ 
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where the essential union is taken with respect to the Lebesgue measure in K+ . 
The bottom of the spectrum of 7i^2'^^'^ is given by 

^^m■.= inf / ^''i, - I,; ■ (5.10) 
Lemma 15 Function (jS.lOp is continuous o?iR+. 

Proof of Lemma Let us fix *p e R+ and e > 0. We will prove that 
+ p) — < e if IpI is small enough. Choose ^ e J)^ such that 



IV^II 



2 



< |. (5.11) 



Let 

We have 

N2 

</> - V = (A|+>' - A|,)^ = E n ^r^A^^^ - A,?j n a^.V'. (5.12) 

/c— 1 i<k j>k 

Let be the unitary Fourier transform in L2{^^^^'^^^^ x S'^,C^"^) defined by 
{TC){^^,^.l, ■ ■ ■ ,qAr2-i) 

W2-1 

:= (2^)3(i-A^^)/2 / ^(a;,yi,...,y^,_i)e"^ '""■'''"dyi • • • dy^,_i. 



K3(N2-1) 

We can rewrite (|5.12p as 

^ _ ^ ^ ^-1 ^ J] A^+P(A,Y+P - A^J n K^^^ (5-13) 

fc— 1 ■i<A; 

where A^, n E Z2 are the operators of muhiphcation by the symbols 

A^ :=^ + "rj?:;^^"r , (5.14) 
Ar2-i 

^ ^ (5.15) 
Pn. := ^^*P + qfe-i, k = 2,...,N2. 
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The matrix-functions (j5.14p are uniformly continuous in *p. Thus by (j5.13p 

N2 

ll</'-'/'II^V2(K3<«.-i,,s.,c^«.) s^cg||A^+P-A:fjl^v2^,/i/= —^ 0. (5.16) 



We write 



(5.17) 



The second and third terms on the r. h. s. of (|5.17p tend to zero as |p| ^ 
according to (|5.16p and (|2.8p . The last term also tends to zero for small |p|, 
since the symbol of the difference is 



From ([Sre]) and (ISTT]) follows that 



■,2 



|, (5.18) 



11^112 11011^ 
if IpI is small enough. Hence by (|5.1ip and (|5.18p for any e > 

|m(^5 + P)-/^(*P)| <e 

for IpI small enough. • 

Now we prove that /i is semibounded from below and tends to infinity as — > 
00. This, together with (15. 9p and Lemma [151 implies that the spectrum of 
7^22'^^ is purely essential and is concentrated on a semi-axis. Proposition 1141 
will be thus proved. 

According to (jT^ for j = 2 and (|2.10p we have 

(H^f^ij,^) > C,{J2 v/-A„ + ™2^,V)-C2||^||^ 

"^^2 (5.19) 
for any ^ e P° ® AnH^^^{R^ ,C^). 

neZ2 

Since all the operators corresponding to the quadratic forms involved in (|5.19p 
commute with J^q^^J^q, it follows from (15. 8p that for almost all ^ 

(H^-^-'PV,^)* ^ Ci{ J2 ^Jvl^TalTi,,Ti,), - C^ml (5.20) 

neZ2 

holds for every G !D^, where p„ are defined in (|5.15p . Thus \i is semibounded 
from below. Since by (|5.15p 

ne22 
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there exists n ^ Z2 such that 




and hence 



neZ2 



1. 

N2 



Thus the r. h. s. of (|5.20p tends to infinity as *p 



00. • 



6 Absence of Gaps 




(6.1) 



Let us first fix a decomposition Z on which the minimum is attained in (|2.20p . 
Following the general strategy of [12], we will prove that for any irreducible 
representations {Di, Ei; D2, E2) ~< {D,E) any 



belongs to cr{Ti.j^' )■ This will imply (|6.ip according to the definition (|2.19p . 
Let 



We will use the notation and results of Section \5\ The following lemma is a 
slight modification of Theorem 8.11 of [12] and is proved along the same lines: 

Lemma 16 Let A be a selfadjoint operator in a Hilbert space ^) and U{'^) he 
a continuous representation of a compact group T by unitary operators in 9j 
such that J7(7)Dom>l C DomA and U{j)A = AU^-y) for any 7 e F. Then 
for any irreducible (matrix) representation D of T the corresponding subspace 
reduces A. For every A € a{A^) where A^ is the reduced operator and 
every e > there exists a D-generating subspace G o/DomA such that 



Remark 1 7 Recall that a subspace G of 9j is called D-generating if the op- 
erator U{'j)\G is unitary in G for all ^ € T and there exists an orthonormal 
base in G such that for every 7 G F the operator U{'y)\G is represented by the 
matrix D{'-f). 

Proof of Lemma 1161 Let r be the dimension of the representation D : j 
{Dik{'j))i f._j^- Let us introduce in S) the bounded operators Pik by 



z 



xi{Z,Di,Ei) + H2{Z,D2,E2) 



Ai := X- Ki{Z,Di,Ei) ^ K2{Z,D2,E2). 



(6.2) 



\\Au — Xu\\ ^ £||it||; for all u E G. 



Pik :=r^Afc(7)C/(7)rfM(7), l,k = l,...,r, 
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where ^ is the invariant probabihty measure on F. It is shown in the proof of 
Theorem 8.11 of jT2] that Pu are orthogonal projections onto mutually orthog- 
onal subspaces of and that 



P^^^Pzi. (6.3) 



1=1 

In fact, Pu is the projection on the subspace of function which belong to the 
row of the representation D. Moreover, Pik is a partial isometry between 
Pkk^ and Piif). Since A G (t{A^), there exists a vector uq G Dom^^ such that 

\\A^uq ~ XuqW sC e||uo||- 

It follows from (|6.3p that there exists I G {1, . . . ,r} such that ||PnMo|| ^ . 
We can thus define ui := PuUo/\\PuUo\\ and then Uk ■= Pkiui for k = 1, . . . , r. 
The subspace G spanned by {uk}l.^i satisfies the statement of the lemma. • 
Let 

rj := dim(i:)j (g> Ej), j = 1, 2. (6.4) 

Since >ci{Z, Di, Ei) belongs to the spectrum of (see definition (|2.18p '). 

by Lemma [TBI we can choose a sequence of (Z?i (E) i?i)-generating subspaces 
{GJ^i of Dom(H2 V^') such that for all g G N 

\\H^^f'^g - MZ, Di, Ei)(f>g\\^^ ^ ^q-'U,h,,,, for all G G,. (6.5) 

Analogously, for any *p ^ we can find a sequence {G'f}^i of {D2 E2)- 
generating subspaces of DomTi^j'^^'^ such that 

llH^-^-^'Vf - m(??X|L ^ for all i^'f G G'^. (6.6) 

Moreover, we can choose a set of functions {i^'^iYiLi in DomTi^j'^^'^ with 

in such a way that for every q G N and I = 1, . . . , r2 V'?,/ belongs to the V'^ 
row of the representation (D2 (E" E2) and satisfies (|6.6p . By Proposition [Til 
Lemma |T5l and (|6.2p we can choose *Po hi such a way that 

A*(?5o) - Ai. (6.7) 

We choose Rq > q so that ([^1^ and ([^1^ hold true for aU n,j = 1,...,N, 
n < j with 

e := g^i(7Vi + l)-!^^^-^/'^^/' (C2 + |Ai| + 2)-'^\ (6.8) 
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where Ni^2 are the numbers of particles in Z1.2, and Ci_2 are the constants in 



(|2.9p for i — 2, and so that for some orthonormal base {4iq,kYk=i of Gq 



L2(R3"i,C-*™i) 



4d2 



(6.9) 



where is the dimension of E', ri_2 are defined in (|6.4|) . and the constant 
t'o > depending only on i?, E\ , E2 will be specified later in the proof of 
Lemma [21] 

By Assumption [5] and Lemma [T51 we can choose a sequence of positive numbers 
tending to zero in such a way that 



|/x(q3)-AiU<7-i foraU *Pe[«Po,'P 



+ OqJ, 



1 

2^ 



{^0+5qf5qCR^<q-\ 



(6.10) 
(6.11) 



where Cr^ is the constant in ()2.1ip . and 

^(^5o + ^,)%.l.i?3<^^ 
Let us choose a function fq e _L2(M+) with supp/g C [*Po,?5o + 5q] so that 

¥0 



|/,(*P)|Vrfqi = i. 



(6.12) 



(6.13) 



Let 



V'9,;(yo, • • ■ ,yw2-i) 
1 



(2^)^ 



/,('P)^f,(^,yi, . . . ,y^,_i)*p2d^dq3, 



(6.14) 



where {yo, . . . , yAr2-i} and {x„} 

are related by (|5.ip and (|5.2p . It follows 
from l|6.13p and the choice of tlf^^ that 



V'9,i||fl2 



1, l^l. 



(6.15) 



and that ■(/'g,; belongs to the I row of {D2 <^ £2). Clearly the linear subspace 
Gq spanned by {4>q,iYi^i is a {D2 ® i?2)-generating subspace of DomH^j'^^- 

Lemma 18 For any n ^ Z2 and ^ e Gq with WiIjW — 1 the one-particle density 



satisfies 



) (dx„j • • -dXn )/dx„ 
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Proof. By (OI)) 

l|P'0,n|Uoo(R3) 5^ (27r)"^/^||pV',ri||Li(K3) 



(27r)s 



g-ip(yo+r„)g-»q3^yoj^(q3) 



(6-16) 



X duj duj dyo rfyi • • • dypj^^i 



dp, 



where r„ x„ — yo, see (|5.2p . Integrating the r. h. s. of (|6.16p in yo wc obtain 
(27r)^5(p + *Ptj — from all the factors involving yo. Estimating the absolute 
value of the integral by the integral of absolute value and taking into account 
that / S{p + . . . )dp = 1 we get 

\\p^Alu^^)^t7z-^,\ / / / / \!,m\\um\ 

X |^f(c.,yi,...,yjv,-i)||vf('^,yi,...,yjv.-i)|r$' (6-17) 

xd^id<^dujd^dyi---dyN,-i ^ -^MVo + Sg)^Sg, 

[Zn 

where at the last step we have used Schwarz inequality and WtpW = 1. The 
formal calculation (|6.16p — (|6.17p is justified by the fact that the integral over 
M.^^^ can be considered as a limit of integrals over expanding finite volumes, 
since V e L2(K^^')- • 

Corollary 19 For any W G L2(M^), n e Z2, and ij: € Gq with ^ 1 we 
have 



.3JV2 



2 1 

|VF(X„)?A(X„1, . . . ,X„„J| dx„, •••dx„„^ ^ TT^iVo + Sq)'^6q\\W\\ 



Let Fq be a subspace of at spanned by the functions 

ipq,k.l{xi, ■ ■ ■ ,XAr) := 0,,fc(Xj : j e Zi) (g) ■lpq,l{Xn : U £ Z2), 

(6.18) 

fc=l,...,ri, ; = l,...,r2, 

where {'/>g,fe}^Li ^'^^ {ipq.i}^ti orthonormal bases of Gq and Gq, respec- 
tively. We obviously have \\(pq,kA\\ l^(^^3n ^c-^i^ = 1- 

Lemma 20 For any g e N C DomHAr. For any if G Fq 

\UHn-XU\ <5g-V;/'r^/'||(p||. 
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Proof. It is enough to show that the functions (|6.18p belong to DomTijv and 
satisfy 



Indeed, by triangle and Cauchy inequalities for 



(6.19) 



ri r-2 
fc=l 1=1 



(6.20) 



we have 



k,l 



k=l 1=1 



^ SUp||(Hjv - A)(^5,fe,i|kl^^?-2^^l!¥'ll 



The domain of Hn can be characterized as the set of functions ^ from 
the form domain §5 A„i/^/^(IR'^, C^) on which the sesquilinear form 

n— 1 

{Hn^,,') is a bounded linear functional in S^n- Functions (|6.18p belong 
to «) A„iJi/2(]a3 ^4) by construction. By ([23]), and dM]) we have 



n<j 



(6.21) 



The sesquilinear forms {{Hzs + 'Hz,2)'Pq,k,h ■} are bounded linear functional 
over L2(] 



p3Af ,^4' 



since (f)q,k G Dom(7i^Y^^) ^^"^ V^g ' ^ DomiJ^ 



More- 



over, by (|6.5p 



q,k 



and by (EH), (ET]), (QUI) . (ICT) . and (ICT) 



(6.22) 



In view of (|6.2ip — (|6.22p and (|6.2p . to prove that lySg,*;,; e DomHN and that 
(|6.19p holds true it is enough to obtain that 



n<3 



(6.23) 
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To do this, we first note that by (|2.12p . (|2.13p . and Cauchy inequahty 



neZ2 



<e(iVi + l) Y \\\Dn\h,,l\\ ^ + W [ \\\Dn\hq,l 

By mUilf . and ([Ol), 

Y ^ C^'{\\in^f' - Xi)%,i\\ +G2 + |Ai|) 

neZ2 

<Cfi(C2 + |Ai| + 2q-i). 
Thus by ([621, dOS]) and for g ^ 1 we obtain 



(6.24) 



(6.25) 



Y ^n-^{|x.|>i?<,} + Y t^«J^{|x„-x,|>H, 



- n£Z2 

Now the scalar functions 



n<j 
n#3 



} 'Pq,kd 



(6.26) 



V„.q{-x.) := |t4(x)|/{|x|^7?,}(x) and C/„j.g(x) := |C/„j(x)|/{|x|^flj(x) 

(6.27) 

are square integrable by (|2.1ip . By Corollary [HI for n G Z2 



(6.28) 



and for n < j, n^j 



2 1 

\\Unj,q'Pq,kd\\^ ^ SUp 1 1 C/„j- g (• -z)?/-,.;]! S$ ^(^g («Po + ^g)^ 1 1 1 1 1, (R3) . 



Hence by (ISTTTD . dOSl) . (PTTjl . and (jCTj) 



(6.29) 



(6.30) 



It remains to add (|6.26p and (|6.30p to obtain (|6.23p . finishing the proof of the 
lemma. • 

The subspace Fg spanned by the functions (|6.18p is Di ® Ei ® D2 ® E2- 
generating. Since it is a sum of Dx ® I?2-generating subspaces, and 
(£>!, Ei;D2, E2) -< {D,E), Fq contains some nontrivial D-generating subspace. 

Hence the subspace Kq := P^Fq is not equal to {0} and is contained in Fq. 
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Lemma 2 1 There exists a constant Ce > such that for every q 

IIP^^II ^ CbII^II, for all ipeFg. (6.31) 
Proof. The projector (|2.15p can be written as 



^"" = 1^ ^ E ^^^W'^- (6.32) 

"TTGnfxnf " 7rGn\(nf xnf ) 



We will denote the first term in (|6.32p by , and the second by . Then 

llP^^f = P^^) = Q^^) + R^ip). (6.33) 
Relation {Di, Ei, D2, E2) -< {D,E) implies that the representation E\I{f x Ilf 
is unitarily equivalent to a sum © niE^'^\ where > are muhiplicities of the 

i=0 

irreducible representations of the group nf X nf with £;(o) = Ei(g>E2. 
For the corresponding characters this gives 



^£;(7r) = E^^^^'-'W' foraU vrenfxnf. 

i=0 

Hence 

fe 

1=0 

where i^i > and Pi is the projector corresponding to the representation E'^^^ 
By construction, Pt^Lp = Lp for any ip € i^g, hence Piip = for i = 1, . . . , fc. Thus 
for any if € Fq 

{p,Q''p) ^uoyf, iyo>0. (6.34) 

We will now estimate the second term on the r. h. s. of ()6.33p . For any n E Z2 
and any ip E Gq with \\t/j\\ = 1 by Corollary [T9l and (|6.12p we have 

For any functions (|6.18p and any tt G 11 inequality (16. 9p implies that 



Now if TT e n \ (nf x nf ), then there exists jo G Zi such that ttjo G ^2- Then 
by 



1^0 



(IT h\^,\<R,}\'P'l,k,l\,Tr\Vq^lj\) ^ (l¥'g,fe,iU{lx,J<i?^,}7r|(^^_jj|) ^ . 
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Thus 

\{^,.k,i,7rip^rkj)\^:r^. ^en\(nfxnf). (6.36) 

Any ip £ Fq can be written as (|6.20p . By (|6.36p and Cauchy inequality for any 

7ren\(nf xnf) 

|(<P,7r(^)| =^ \^ki\\c]^\\{Vq,k,hTr^g^],j)\ ^ ^ll^ll^- (6-37) 

k,l,kj ^ 

Since the number of elements of n\ (Ilf x Ilf ) does not exceed du and for any 
TT |^_E(7r)| < as a trace of unitary matrix of dimension ds, (|6.37[) implies 
that 

By (|03l) and (fOi]) wc conclude that (lOTj) holds with Cb = ^/uo/2. • 
Lemmata [20] and [21] imply that ig := P^Kq is a nontrivial subspace of 
DomTi^'^ and for every / = P^ip e Lq 

||(7^^^^-A)/|| ||(7^jv-A)^|| ^5q-\lrlM^5q-\lrlcE'\\fl qeE. 
This implies that A G (t{H^'^), and thus finishes the proof of Thcorem[Sl 
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